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ABSTRACT 

We use the Expanded Very Large Array to image radio continuum emission from local luminous 
and ultraluminous infrared galaxies (LIRGs and ULIRGs) in 1 GHz windows centered at 4.7, 7.2, 29, 
and 36 GHz. This allows us to probe the integrated radio spectral energy distribution (SED) of the 
most energetic galaxies in the local universe. The 4-8 GHz flux densities agree well with previous 
measurements. They yield spectral indices a ~ —0.67 (where F^, oc v") with ±0.15 (Icr) scatter, typical 
of nonthermal (synchrotron) emission from star-forming galaxies. The contrast of our 4-8 GHz data 
with literature 1.5 and 8.4 GHz flux densities gives further evidence for curvature of the radio SED 
of U/LIRGs. The SED appears flatter near ~ 1 GHz than near ~ 6 GHz, suggesting significant 
optical depth effects at the lower frequencies. The high frequency (28-37 GHz) flux densities are low 
compared to extrapolations from the 4-8 GHz data. We confirm and extend to higher frequency a 
previously observed deficit of high frequency radio emission for luminous starburst galaxies. 



1. INTRODUCTION 

Most very luminous galaxies at z ~ emit the bulk of 
their luminosity in the infrared (IR) . These luminous and 
ultraluminous infrared galaxies (LIRGs and ULIRGs, in- 
frared luminosity Ljr [8 — lOOO/um] > 10"'^^ and Lir > 
10^^ Lq) make up the majority of galaxies with Lboi ^ 
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These systems often correspond to major 



mergers and compared to normal star-forming disk galax- 
ies, U/LIRGs exhibit enhanced, concentrated molecular 
gas content, high star formation efficiency, increased ac- 
tive galact ic nuclei (AGN) activity, and disturbed mor- 
phologies (jSanders fc Mirabel 119961) . Even the closest 
ULIRGs are > 100 Mpc away with compact radi o, IR , 
and millimeter line emission ("e.g.. iLo nsdalc ct al.l 119931 : 
iDownes fc Solomonlll998t iSakamotoet al...2008. ). so that 
high angular resolution (such as that obtained by the 
EVLA) is required to study their structure. 

Nonthermal (synchrotron) radio emission is well- 
established t o track infrared e mission (the "radio-FIR 
corre lation," iHelou et al.lll985l : lCondonlll99l lYun et al.l 
l2001h while thermal (bremsstrahlung) radio emission at 
v ^ i GHz gives an unobscured view of ionizing photon 
production. These emission mechanisms, the indifference 
of radio waves to dust, and the resolving power of cm- 
wave interferometers make the radio continuum a natural 
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tool to dissect the detailed structure of U/LIRGs. De- 
viations from the radio-IR correlation, radio morphol- 
ogy, and brightness can all help distinguish the pres- 
ence and contribution of AGN to the overall luminos- 
ity, though not necess arily uniquely ([Condon et al. 199^ 
'Lons dale erallll995t fSaiina et al.l [20081: iLonsdale et atl 
2006, and references in the latter). 

This letter reports first results from our Expanded Very 
Large Array (EVLA) program to explore the structure of 
star formation in U/LIRGs using radio continuum emis- 
sion. We are imaging 22 of the most luminous nearby 
U/LIRGs (LiR > 10"-6 Lq) in the C (4-8 GHz) band and 
Ka (26.5-40 GHz) bands. Here we report observations 
from the EVLA's C and D configurations. The EVLA's 
continuum sensitivity makes it possible to achieve good 
signal-to-noise ratio, even at high frequencies, with mod- 
est time on source. The large instantaneous bandwidth 
allows us to probe the radio spectral index, a (defined so 
that Fi, ex v"), both within and between bands. 

Here explore variations in the radio spectral energy 
distribution (SED) of our targets and in the radio-IR 
correlation. In a simple picture dominated by star for- 
mation we would expect the C band to be dominated 
by synchrotron emission with spectral index a ~ —0.8 
while the Ka contains an approximately equal mixture 
of nonthermal emission and thermal emission, which has 
I' = —0.1. For systems dominated by an AGN the spec- 
tral index may show wider variation depending on the 
orientation of the AGN or injection spectrum of electron 
energies. We close by showing overlays of the Ka-band 
emission on HST imaging. These highlight the EVLA's 
improved ability to trace embedded star formation at 
high resolution. 

2. OBSERVATIONS AND DATA REDUCTION 

Our EVLA resident shared risk observing program, 
AL746, is imaging radio continuum emission from nearby 
LIRGs and ULIRGs. The sample, listed in Table [l] con- 
sists of 22 of the most luminous members of the IRAS Re- 
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vised Bright Galaxy Sample (RBGS [Sanders et al.|[2003l) 
northern enough to be observed by the EVLA. Distances 
span 50-300 Mpc, with a median of 177 Mpc, and IR 
luminosities go from log^Q Lir [L0] = 11.6-12.(0. 

The whole program consists of imaging each target at 
four frequencies, v ~ 4.7, 7.2, 29, and 36 GHz, with 
a bandwidth of 1 GHz centered at each frequency. We 
observe the two low frequency windows (together the "C- 
band" data) and two high frequencies (the "Ka-band" 
data) simultaneously. Here we report the results from 
the two most compact EVLA configurations (C and D). 
We observed each target for ^ 5 minutes in each band 
and configuration. 

We reduced the data using the Common Astronomy 
Software Applications (CASA) package developed by 
NRAO. We adopted standard approaches to flag bad 
data (radio frequency interference is a particular issue in 
the C band) and calibrate the bandpass, phase, and am- 
plitude response of individual antennas. We iteratively 
self-calibrated the data whenever possible. 

We combined all data for each band into a single image 
using CASA's multi-frequency synthesis (MPS) capabili- 
ties (see Sault & Wieringa 1994). This image reconstruc- 
tion mode synthesizes interferometer observations across 
a wide bandwidth. It accounts for changes in the syn- 
thesized beam as a function of frequency and produces 
images of intensity for a fiducial frequency and maps of 
the spectral index. For our data, these fiducial frequen- 
cies are 5.95 GHz for the C band and 32.5 GHz for the 
Ka band. Note that we do not actually observe at either 
of these frequencies, they are merely a intermediate be- 
tween our two separated 1 GHz windows. At C-band, we 
deconvolved (CLEAN) the image using components that 
have both a flux density and a spectral index. At Ka 
band we use components that have only flux density at 
the fiducial frequency because we expect a flatter spec- 
tral index, have a lower signal-to-noise ratio, and have a 
smaller fractional bandwidth. 

After reduction, we have maps of intensity for C and 
Ka band and spectral index at C band. Synthesized 
beams are ~ 8" at C band and ~ 1" at Ka band. The 
noise in a C-band image is usually k, 20 /iJy beam~^, 
though for our brightest targets calibration issues some- 
times limit us to about twice this value. A typical noise 
in a Ka-band image is w 50 /iJy beam"-'^. All targets are 
detected, with peak SNR from 200 to 6.3 x 10^ at C band 
(median 670) and 12 to 500 at Ka band (median 50). 

We measure C-band flux densities in the following way: 
we estimate the RMS noise away from the source. We 
then identify the region of contiguous SNR > 10 emission 
that contains the galaxy. We integrate over this area to 
derive the flux density. We take the intensity-weighted 
mean of the spectral index (ac) = J F^a/ J F^, derived 
by the MFS CLEAN as the spectral index for the source. 
We derive Ka-band flux densities by integrating over the 
same area used to derive the C-band flux densities. Over 
this matched area, two galaxies have integrated SNR < 3 
at Ka and we omit them from all Ka-band related anal- 
ysis. 



® This refers to estimates of the "total" (8- lOO Oum) infrared 
luminosity eonstructed from IRAS flux densities bv lSanders et aP 
p003 ) following Sanders fc Mirabel (1996 ). 



3. RESULTS 
3.L A Typical Synchrotron Spectrum at C'-Band 

The left panel in Figure [1] compares our measured C- 
band flux densities at 5.95 GHz (x-axis) to previous mea- 
surements at 8.4 GHz (black cir cles) and 1.5 G Hz (open 
squares) both drawn from lCondon et al.l ()1991l ). Table [T] 
reports internal C-band spectral indices, ac- Our median 
ac is -0.67 with la scatter « 0.15. The SNR of the data 
are very high but an uncertainty of even a few percent in 
the relative calibrations of the two parts of the C-band 
(e.g., due to calibration problems, radio frequency inter- 
ference, uncertainties in the primary calibrator model, or 
CLEAN issues) will lead to uncertainties '^0.1. 

We expect nonthermal (synchrotron) emission to dom- 
inate the C-band. Our measured a — —0.67 is close to 
a normal synchrotron spectral index for a star-forming 
galaxy, a '« -0.8 (Condon 1992) though AGN can 
produce similar indices. For ac = —0.67, Figure [1] 
shows that our C-band flux densities agree with the X- 
band measurements of iCondon et al.l (|1991f) well. The 
median ratio after correcting to the same band using 
ac — —0.67 is « 1.02 and the scatter is 14%, indicating 
very good agreement. Thus our C-band measurements 
from « 4.2-5.2 and 6.7-7.7 GHz and previous measure- 
ments at 8.44 GHz all exhibit good consistency with a 
typical synchrotron spectrum. 

3.2. Departures From a Simple SED at Low and High v 

In the simple two-component model outlined in the 
introduction synchrotron emission with a ~ —0.8 dom- 
inates below V ^ Vd GHz. Above this frequency, the 
radio continuum results from a mixture of synchrotron 
and thermal emission, which has a = —0.1. In a nor- 
mal star-forming galaxy the ratio of these two processes 
may be approximately fixed. Figure [2] shows the radio 
SEDs of our targets, which exhibit significant departures 
from this picture at both high and low frequency. We 
normalize eac h SED at 5.95 GHz, show literatur e flux 
densit ies from ICondon et al.l ()1991l ) and Clemens et al.l 
(J2608f l. and over plot the s i mple t wo-component, optically 
thin model from iCondonI ()1992[ ). The lines through our 
C-band points show our measured ac- 

At low frequencies, the 1.5 GHz (L band) flux densi- 
ties lie mostly below an extrapolation of the synchrotron 
emission measured in C-band, which appears as a dashed 
line in Figure [2l This can be seen from the left panel in 
Figure [1] where the 1.5 GHz points (open squares) tend 
to scatter below the value expected for our measured ac 
(dashed). Our C-band flux densities and ac demonstrate 
normal synchrotron emission at z^ > 4 GHz. The low L- 
band flux densities imply that the radio SED becomes 
flatter (a approaches 0) near v ^\ GHz. This effect was 
previously inferred from unusually flat (near 0) spectra l 
indices contrasting 1.5-to-8.4 GHz (jCondon et al.in"991f ) 
an d directly observed v ia long wavelength observations 
bv lClemenseFaLl (|2010t ). 

For a normal star-forming galaxy, we expect the radio 
SED to contain an increasing fraction of thermal emis- 
sion (a = —0.1) moving toward higher frequencies. By 
~ 30 GHz (the Ka band) , the SED of such a galaxy would 
contain roughly equal parts thermal and non-thermal 
emission with flatter a than at 5.95 GHz. Figure [2] shows 
that this simple expectation does not hold for our data. 
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Fig. 1. — (left) Flux densities from [Uondon et ai.l 119911) (j/-axis) as a function of our 5.95 GHz (C-band) flux density (x-axis). FiUed 
circles show X-band (8.4 GHz) flux densities. Open squares show L-band (1.5 GHz) flux densities. Dashed lines show the expected flux 
density given our C-band flux density and a spectral index a = —0.67. Our C-band flux densities agree well with the X-band values. The 
L-band flux densities tend to lie below the expected value, indicating curvature in the radio SED in the range 1-5 GHz. (right) Ratio of 
far-infrared to radio emission (q) for Ka (open diamonds), C (black circles), and L (open squares) bands. All three bands correlate with 
IR emission, as expected, with the lowest scatter in the C-band. The highest luminosity source, with an inverted spectrum, is UGC 08058 
(Markarian 231), which is dominated by an AGN. Red outlines highlight source and two other outliers discussed in the text. 



Our Ka-band flux densities tend to lie below the sim- 
ple model: the model predicts a typical C-to-Ka ratio 
Ri 2.4, while our median ratio is ~ 3.f5 with la scatter 
0.75. The radio spectrum is significantly steeper than 
expected moving to higher frequencies. 

A steep SED at high f requencies has been p reviously 
observed in U/LlRGs bv lClemens et a l.' (2008^, and the 
dwarf starburst NGC 1569 (iLisenfeTd et al. 2004). If 
there is a significant thermal contribution at i/ '^ 30 GHz, 
and there must be based on the high level of star for- 
mation in these systems, then the synchrotron spectral 
index must be steepening. iCleme ns et al.' ('2008') discuss 
possible mechanisms for this and Cl e mens ct al. (201(3) 
revisit the topic, while iLisenfeld et aLl ()2004f ) discuss pos- 
sible explanations for a similar observation in the star- 
burst NGC 1569. Allowing that the thermal emission 
should be present, the synchrotron emission must be de- 
pressed at high frequency. Depressed synchrotron emis- 
sion could result from synchrotron and inverse Compton 
losses in aging starbursts or regions with high magnetic 
fields. However, there is no evidence that we observe all 
of these galaxies at a preferential time and these expla- 
nations appear inadequate to match the magnitude of 
the effect. Alternatively, the injection spectrum of elec- 
tron energies could differ from those in normal galaxies. 
This is possible for an AGN, but less plausible for ra- 
dio continuum emission originating from star formation. 
We suspect that a modified injection spectrum leads to a 
deficit of high frequency synchrotron emission and drives 
the high-frequency flux deflcit. This is a topic that we 
will investigate further with our full high resolution data. 

Depressed synchrotron at high frequencies may be a 
way to distinguish systems where the radio is dominated 
by AGN. We expect a minimum level of thermal emission 
in systems with radio emission due mostly to star forma- 
tion, so AGN likely provide most of the radio emission 



in the steepest-spectrum sources. 

3.3. Radio-Infrared Correlation 

Despite the surprises in the radio SED, the radio and 
infrared fluxes still track one another well. We calcu- 
late th e far-infrared-to-radio ratio deflned by Helou et al] 
<|1985f l as q^ = logio (i^FiR/5'i.) with i^piR the far- 
IR flujO in units of 3.75 x lO^^ W m'^ and S^ in 
W m^-^ Hz^^. The median 55.95 = 2.8 with a scatter 
of sa 0.16 dex. The Ka band shows comparable scatter. 
The L-band (/1.5 has ^ 50% more scatter, '^ 0.25 dex. 
lYun et ahl (|200lr ) observed 0.26 dex scatter in the L-band 
radio-IR correlation and 0.33 dex scatter for LIRGs and 
ULIRGs. The smaller overall scatter in our sample may 
simply be a result of small number statistics. We inter- 
pret the lower scatter at high frequencies in our sample 
as a manifestation of the same phenomenon that leads 
to higher scatter in (71.5 at high luminosities. The curved 
SED at low frequencies in luminous IR galaxies adds scat- 
ter to the radio-IR correlation as physical processes other 
than just star formation become important. At higher 
frequencies and in lower luminosity galaxies the effects 
are less important. 

The correlation among the radio bands is stronger than 
the correlation with the IR, with the strongest correlation 
among the two synchrotron dominated bands: 55.95 and 
932.5 rank correlate with 0.67; (75.95 and gi.5 exhibit a 
rank correlation coefficient of 0.81; 91.5 and (732.5 have a 
rank correlation 0.57. 

Three galaxies lie notably off the radio-IR correlation. 
The radio emission from UGC 08058 (Markarian 231) is 
much (4.5(7) brighter than one would derive from the IR 
emission. UGC 05101 shows a weaker (2.5(j) excess in the 

^'^ To calculate q we use the IRAS-defined far-IR flux, which 
captures flux from 40-120/xm, not the integrated IR flux quoted 
elsewhere in this paper. 
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Fig. 2. — Normalized spectral energy distributions of our targets, sorted by IR luminosity from top left to bottom right. Gray regions show 
our observed frequencies and circles indicate our measured fiux densities. Lines show ac across our C-band measurement. We normalize 
all values to the 5.95 GHz value and plot literature measurements as open symbols. The dashed line shows a simple two-component model 
mixin g optically thin synchrotron (a = —0.8) with thermal emission (o = —0.1) at levels appropriate for a star-forming galaxy (Condo^ 
119921) . Our C-band measurements show an approximately normal synchrotron shape in most sources, but U/LIRGs depart from this simple 
model at low and high frequencies. Our new Ka measurements confirm a systematic deficit of high-frequency emission, extending this 
result to higher frequencies [u ~ 36 GHz). 

same direction. Bot h galaxies show evidence fo r AGN in 
VLBI observations fLonsd ale et al.]|1993ll2003l) . Though 
the contrib ution of the AGN to t he bofometric luminosity 
is debated, iVeilleux et alj (|2009D estimate the fraction of 
total IR emission due to the AGN to be ~ 70 ± 15% in 
UGC 08058 and ~ 35 ±15% in UGC 05101 IArmus"^ran 
(see [20071 for an even lower estimate). 

IRASF 08572-^3915 shows unusually (2.5cr) strong IR 
emission. This is our faintest target a nd its IR-excess 
has been noted before (see discussion in lYun et al.ll2004() 
with ~ 25% of the IR emission in this galaxy asso- 
ciated with a very hot T > 300 K dust component 
(jArmus et al.l 120071 the largest such fraction in their 
sample). The flat spectral index at low frequencies 
and the high IR-to-radio ratio in this galaxy could be 
naively interpreted to indicate a very young starburst 



(e.g. IRoussel et alll2006D . For example, this galaxy sits 
next to the yo ungest-age model point in the diagrams of 
iBressan et al.l (f2002.') . However, based on near- and mid- 
IR spectroscopv (lArmus et al.ll2007t i Risaliti et al.ll2010D 
and X-Ray hardness ([Iwasawa et al.ll2011[ ) much of the 
luminosity is thought to come from a deeply embedded 
AGN in the northwest nucleus. In either case, a lumi- 
nous, compact, deeply embedded object drives the lu- 
minosity and results unusually strong IR emission. We 
detect it at Ka band, but not at > 3a over the C-band 
aperture so that further study will have to wait for the 
the better SNR and smaller matched apertures from our 
full data set. 



3.4. Tracing Embedded Star Formation 
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Fig. 3.— EVLA Ka-band contours superimposed on grey scale HST ACS/WFC I-band (F814W) images of CGCG 448-020 (II Zw 096, 
left) and NGC 3690 (right), both on a logarithmic stretch. The contours are displayed in six evenly spaced log steps between 10-100% of 
peak intensity for II Zw 096 and 2-100% of the peak intensity for NGC 3690. 



Figure |3] shows our Ka-band emission superimposed 
on /-band HST images for two targets, CGCG 448-020 
(II Zw 096, left) and NGC 3690 (right). The images show 
the utihty of combining high-resohition radio data with 
the muhiwavelength data assembled by t he Great Ob- 
serva tory All-sky LIRG Survey (GOALS, lArmus et all 
|2009() to assess the distribution and strength of embed- 
ded star formation and AGN activity. In II Zw 096, the 
brightest radio continuum point does not have an ob- 
vious optical counterpart. This bright radio emission 
coincides with the peak of t he IR emission measured 
by Spitzer (jlnami et al.l 120101 ) . underscoring that this is 
a star-forming peak without an obvious optical signa- 
ture. In NGC 3960 (right), several of the radio sources 
have optical counterparts, but without the radio emis- 
sion it is difficult to determine which optical knots har- 
bor luminous embedded activity. Here too, the 32.5 GHz 
continuum peaks coincide with the bright IR emission 
(jAlonso-Herrero et al.ll20io() but with higher angular res- 
olution. 

3.5. Summary 

We present new measurements of radio continuum 
emission from 22 U/LIRGs across 1 GHz windows cen- 
tered at 4.7, 7.2, 29, and 36 GHz. Our observations 
are consistent with typical synchrotron emission at 4.2- 
7.7 GHz. Comparing these normal C-band results with 
new high-frequency 29-36 GHz fluxes and literature mea- 
surements at L, X, and K-bands, we see that the radio 



SEDs of U/LIRGs systematically deviate from a simple 
two-component model at both low and high frequencies. 
The new EVLA data here provide: 1) observations fill- 
ing in previously unsampled parts of the radio SED; 2) 
a clear observation of "normal" synchrotron behavior in 
the range 4-8 GHz, reinforcing that the radio SED must 
be curving away from simple expectations; and 3) an ob- 
servation of systematically low high-frequency flux den- 
sity at higher frequencies than have been previously stud- 
ied. The radio fiuxes correlate with IR emission as ex- 
pected, with some suggestion that the C- and Ka-bands 
exhibit a more direct link to IR emission than the L- 
band. 

An ongoing extension of this project to all EVLA con- 
figurations will dramatically improve our resolution and 
signal-to-noise ratio in both frequency bands. Combin- 
ing these data with the multiwavelength data available 
for our targets, we expect to test if the star forma- 
tion surface density exhibits any maximum value (e.g., 
the limit set by mo mentum-driven winds suggested by 
iMurrav et al.l 120051 ) , derive characteristic sizes for the 
star-forming disks, map the spectral index at < 0.5" 
resolution, and explore the implications of the bright- 
ness temperature and radio SED for the dominant radio 
emission mechanism. 
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Radio Spectrum of U/LIRGS 



TABLE 1 
Integrated Properties 



Galaxy 



Dist.'' LiR^' 

[Mpc] [logio Lq\ 
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[mjy] 



ac 



S'32.5'' 

[mJy] 



5i.49 (C91)d 
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S8.44 (C91)d 
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Si5 (C08)d 
[mJy] 



CGCG 436-030 
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11.69 


18.6 + 0.04 
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43.1 


12.7 


8.73 


CGCG 448-020 (II Zw 
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14.6 + 0.06 
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25.4 + 0.05 
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7.3 + 0.6 
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9.5 + 0.04 
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IRAS 21101-1-5810 
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9.8 + 0.04 


-0.67 


4.6 + 0.6 










IRAS F01365-1042 
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11.85 


10.0 + 0.04 
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4.1 + 0.5 


17.0 
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IRAS 0857+39 
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4.44 + 0.04 


-0.29 
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4 
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IRAS 1525+36 
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MCG +07-23-019 (Arp 148) 
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^ Hubble flow distance and IR luminosity adopted from 


ISanders et al.l II2003I) foUowine 1 Armus 


et al.| P009|). 









^ Flux density integrated over significant C-band emission. See text. Uncertainties are only statistical, bootstrapped from the map. 

'^ Intensity-weighted in-band spectral index derived from MFS imaging. Sec text. 

"* Literature flux densities from lCondon et al.l 11199 II ) and lClemens et al.l ([2008). 

■^ Signal-to-noise ratio < 3 integrated over the C-band beam. Still detect at Ka resolution. 



